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ABSTRACT 
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The evolution of a two dimensional, incompressible, rapidly decelerating, time- 
dependent viscous flow about a sharp-edged camber is simulated through the use of the 
discrete vortex model. Vorticity is represented by a distribution of discrete point 
vortices. Each vortex is convected in the velocity field, calculated locally using the Biot- 
Savart law. The roll-up of the vortex sheets, the distribution of velocity and pressure 
on the camber, and the drag force are calculated at suitable time intervals for a 
prescribed time-dependent flow. Experiments are carried out in a vertical water tunnel 
partly to measure the drag force and partly to record on a video tape the evolution of 
the v/ake. The measured and calculated characteristics of the flow, such as the growth 
of the wake and the forces acting on the camber are found to be in good agreement. 
Furthermore, the numerical simulation provided a plausible explanation for the cause 
of parachute collapse, a phenomenon which gave impetus to the numerical and 
physical experiments described herein. The numerical model developed during the 
course of the investigation is applicable to any time-dependent flow about two- 
dimensional cambered plates (circular arcs ). 
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I. INTRODUCTION 



A. SEPARATED FLOWS 

The separated flow about blufT bodies has been almost completely unyielding to 
both analysis and numerical simulation for a number of mathematical reasons and 
fundamental fluid dynamic phenomena. Separation gives rise to the formation of free 
shear layers which roll up into vortex rings or counter-rotating vortices. They, in turn, 
interact with each other, with the counter-sign vorticity generated at the base of the 
body, and with the motion of often unknown separation points. The wake becomes 
unsteady even for a steady ambient flow and the problem of the determination of the 
characteristics of the wake becomes coupled to the conditions prevailing upstream of 
the separation points. Evidently, viscosity modifies radically the inviscid flow, which, 
in this case, cannot serve even as a first approximation to the actual flow. The 
boundary layer equations are not applicable beyond the separation points and are, 
therefore, of limited use in bluff-body flow problems. 

The separated unsteady flow situations involving wake return, as in the case of a 
decelerating or oscillating body, are an order of magnitude more complex than those 
where the vortices continuously move away from the body. The net effect of the wake 
return is twofold. Firstly, the proximity of the large vortices dramatically affects the 
boundary layer, outer flow, pressure distribution, and the generation and survival rate 
of the new vorticity. Secondly, the vortices not only give rise to additional separation 
points and/or additional vortices, but also strongly affect the motion of the primary 
vortices. These effects are further compounded by the diffusion and decay of vortices 
and by the three-dimensional nature of the flow. 

The existing finite difference and finite element methods cannot yet treat the high 
Reynolds number flows with sufficient accuracy for a number of reasons. The finite 
difference schemes require a very fine grid, a turbulence model, and a very large 
computer memory. It seems that the modelling of the turbulent stresses in the wake, 
particularly in time-dependent flows will be the major source of difficulty in all future 
calculations. Whether or not it will ever be practical to apply the finite difference and 
finite element methods to high Reynolds number flows is unknown. The inherent 
difficulties are certainly significant enough to warrant exploring other solution 
methods. 
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Certain separated time-dependent flows may be simulated through the use of the 
discrete vortex model (DVM) (see e.g., Chorin 1973; Sarpkaya 1975). The free shear 
layers which emanate from the sides of the body are represented by an assembly of 
discrete vortices. The strength of the elemental vortices are determined through the use 
of the Kutta condition. The use of a suitable convection scheme enables one to march 
in time and to calculate the evolution of the wake, the velocity and pressure 
distributions, and the lift and drag forces acting on the body. The work described 
herein deals with the application of the DVM to decelerating flow about a two- 
dimensional camber. 

B. FLOW ABOUT A CAMBERED PLATE AND PARACHUTE COLLAPSE 

The determination of the deployment sequence of an axisymmetric porous 
parachute and the unsteady aerodynamic loads acting on it present a very complex 
coupled problem. The development of an analytical or numerical model which takes 
into account the effects of porosity, gaps, and variable opening schemes would allow 
numerical experiments on a large class of parachutes, reduce the number of e.xpensive 
field tests to a few judiciously selected ones, and enable the designer to calculate the 
time history of the fall of the parachute and the strength required to survive the 
aerodynamic loads. However, the development of such a model is hampered by a 
number of difficulties. 

The previous models for parachute loads are based by and large on empirical 
assumptions (see e.g., Heinrich and Saari 1987; Mcwey 1972). They rely on the 
observation that families of parachutes open in a characteristic length and seem to 
have aerodynamic properties that relate well to the projected area of the parachute. 
The apparent mass is assumed to be a function of the projected area only and is not a 
function of the prevailing flow characteristics. The vortex sheet analysis was used by 
Klimas (1977) to derive the acceleration-independent apparent mass coefficient for 
arbitrary-shaped axisymmetric surfaces. Muramoto and Garrard (1984) used a 
continuous-source model to predict the steady-state drag of ribbon parachutes. The 
analyses did not, however, deal with the evolution of the unsteady wake and its 
interaction with the canopy. 

It is in view of the foregoing that a fundamental study of the separated time- 
dependent flow about two-dimensional rigid cambered plates was undertaken. Clearly, 
the flow about a rigid cambered plate is considerably simpler than that about a porous. 



axisymmetric. and flexible parachute and the results, regardless of the degree of their 
agreement with corresponding experiments, may not have direct relevance to the 
practical problem under consideration. But the object of this investigation was the 
understanding of the evolution of the wake under controlled conditions rather than to 
provide a design tool. It is hoped that an investigation of this type will reveal the 
underlying physics of the phenomenon (particularly that of the parachute collapse), 
help to interpret the full-scale results and will provide inspiration for the development 
of more general vortex models with which the dynamics of axisymmetric, porous, and 
flexible parachute canopies can be investigated. 
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II. ANALYSIS 



A. TR.ANS FORMATIONS .4ND THE COMPLEX VELOCITY POTENTIAL 

The calculation of the velocity of any one of the vortices and the force acting on 
the body requires a conformal transformation (in which the camber becomes a circle), 
a complex-velocity potential representing the vonices, their images, and the rwo- 
dimensional irrotational flow around the body, and the use of the generalized Blasius 
theorem. 

The flow in the circle plane may be transformed to that about a cambered plate 
through the use of two successive transformations, one from plane to the plane 
and the ocher from the plane to the z plane. These are given by (see Fig. 2.1) 






Figure 2.1 Circle and physical planes 
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Z = . and ^ + m 



(2.1) 



Connbining the two, one has a direct transformation from the plane to the z 
plane as 



z = i^ + m — 



b^ 

^ + m 



( 2 . 2 ) 



It is easy to show that the camber in the z plane is a circular arc. 

The y-axis in the z plane passes through the tips of the camber. It is 
advantageous to locate the origin of the coordinate axes at the geometric center of the 
camber, i.e.. at the center of the circle part of which represents the camber. This is 
easily accomplished by shifting the origin of the coordinate axes by 



where z^' is the x coordinate of the origin of the circle in the z plane. Thus, one has 
b- 

2 + m + z^ , withz^' = -z^ (2.4) 

which transforms the circle in Fig. 2.1a to the physical plane in Fig. 2. Id. Table 1 
summarizes the relationship between m, z^', the included angle of the camber, b, and 
the radius of the camber. 



TABLE I 

SU.MMARY OF THE PARAMETRIC RELATIONSHIP 



m 




2a 


b 


R= 1/m 


cos60 = .5 


-1 


120 


.866 


2 


COS45 = .707 


0. 


180 


.707 




cos30 = .866 


1/V3 


240 


.5 


2/V3 
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The complex potential function W in the circle plane (see I'ig. 2.1a) which 
tlescribes a uniform Ilow U (assumed to be time-dependent) with a doublet at the 
origin to simulate the cylinder, clockwise-rotating vortices (called q-vortices), Fj^p 
counter-clockwise rotating vortices (called p- vortices), and the images of all the p-and 
q-vortices in the circle plane may be written as 

w= - U(; + ^) + Ln(; - ) - !ToPL„(i; - ^ > 




in which Fj^p and i^j.p represent respectively the strength and location of the k-th p- 
vortex, Fj.^ and the strength and location of the k-th q-vortex, and c the radius of 
the cylinder; an overbar indicates a complex conjugate. The need for the separate 
identification of the p-and q-vortices and for the singling out of one of the vortices in 
each shear layer (namely Fqp and Fq^ i.e., the nascent vortices ) will become apparent 
later. 

B. COMPLEX VELOCITIES OF VORTICES 

The convection of the vortices and the calculation of the forces acting on the 
body require the evaluation of the velocities at the vortex centers. For the velocities in 
the circle plane this reduces to subtracting from Fq. (2.5) the complex potential 
corresponding to the vortex for which the velocity components are to be determined 
and evaluating the derivative of the remaining terms at To determine the 

velocities in the physical plane, however, one has to subtract (iFj./2rt) Ln(z — Zj.) from 
Eq. (2.5) or, in terms of^ , the terms (see e.g., Sarpkaya 1967, 1975) 

ir ir 

4^ Ln (i; - C^) + 4^^ Ln (1 - 4 ) (2.6) 

2n ^ 2n 
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It should be noted that the first term in Eq. (2.6) is the complex function 
corresponding to the k.-th vortex in the ^ plane. The second term appears merely as a 
consequence of the transformation used. It is easy to show that for a p-vortex Eq. 
(2.6) may be reduced to (see e.g., Mostafa 1987) 



^^ 1,0 (~b-)(C,, + m) 

27t [(^^p + m)“ + b^ f 



(2.7) 



A similar e.xpression can be written for the q-vortex. 

C. KUTTA CONDITION 

The fact that the flow separates tangentially with a finite velocity at the edges of 
the plate (Kutta condition) may be expressed by requiring 



dW 

— =0 at = C = — m ±ib 

d'; 



( 2 . 8 ) 



Thus, inserting Eq. (2.5) in Eq. (2.8) one has 



if. 



'>jT C — C c^ 

-Op r - _ 

'’t r 
’Op 



irn 

- — Oq 



1 



2^ 



Oq 



^Oq 



m ;p 1 

vJikp L 

^ 2Jt C -C 



1 



m Y 



1 



k=i 



t -kp r 

'’t r 
•kp 



k7,2:t 



kq 



c 

^kq 



-U(l- — ) =0. 



(2.9) 



Equation (2.9) may be decomposed into two parts as 
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1 



1 



I 



I 



ir, 



Op 



2^ 



Op 










+ (~^0 + *'’ o ) = 0 . 



( 2 . 10 ) 



where the terms containing the strength of the nascent vortices represent the velocity 
induced at the tip of the camber by the nascent vortices and the term in parenthesis the 
velocity at the tip due to all other vortices (and their images), the doublet at the center 
of the circle in the ^ plane and the ambient velocity. 

Equation (2.10) represents two coupled equations for the strengths and positions 
of the nascent vortices. Thus, the solution of the said quantities does, in general, 
require an iteration. However, this iteration may be avoided by noting that the 
velocity induced by a nascent vortex at the opposite tip is very small and certainly 
negligible. Thus, Eq. (2.10) for one of the nascent vortices may be reduced to 



-'£o, 

271 



’Oq 






+ (-u^ + iVo) 



= 0 



( 2 . 11 ) 



A similar expression may be written for the other nascent vortex. The use of the Kutta 
condition, as expressed by Eq. (2.11), will be further explained following the discussion 
of the tip velocity. It suffices to note that all nascent vortices satisfying the Kutta 
condition do not yield either the same tip velocity or the same velocity distribution in 
the neighborhood of the tip. There are, in fact, certain preferred positions for the 
nascent vortices which yield physically realistic velocity distributions near the tips of 
the cambered plate. These nascent vortex positions have been determined by Mostafa 
(19S7). 

D. TIP VELOCITY 

According to the Kutta condition the tangential velocity at the tip is finite. It 
may be determined through the use of I'Hopital's rule. 

The velocity at the tip is given by 



dW _ dW di; 
dz d^ dz 



at z = z„ ± 2ib 

t O 



( 2 . 12 ) 
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The use of rHopital's rule then yields, 



dW _ d^W ib 

dz 6^ ^ 2 ^ 

z = z. 



(2.13) 



Equation (2.13) yields the desired finite tip velocities. 

E. TIME DEPENDENT-FORCES 

The force acting on the body in the physical plane may be calculated either 
through the use of the pressure distribution or through the use of the rate of change of 
impulse. 

Bernoulli's equation for unsteady flow is given by 




2 

il ^ - 



) - ( 



7 V,- 2 

- I — ds = f(t) 
P 2 / 5t 



(2.14) 



where the indices indicate two points on the body in the physical plane. Since there is 
no pressure drop across the shear layer and since the integral term in Eq. (2.14) is zero 
at the tip (i.e., ds = 0), one has 



f(t) = 





(2.15) 



where and represent the tangental velocities on the upstream and downstream 
faces of the tip. It is important to note that f (t) in Eq. (2.15) is also the time rate of 
change of circulation, i.e., the rate at which vorticity is shed into the wake from the tip 
of the cambered plate. 

The normalized form of Bernoulli's equation between any two points m and n 
then becomes 




V 2- V 2 
tl ^2 



U 



• 2 



4- 



V 2 - V 2 
n ni 



u 



2 

o 



d " 




(2.16) 
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The integration of the differential pressure between the upstream and downstream faces 
of the camber yields the force components in the x and y directions, i.e., the drag and 
lift forces. 

The force acting on the body can also be calculated through the rate of change of 
impulse. It is given by 



F = Anpc^U (1 - 



m" d 



ny 



(2.17) 



which may be written as 



+ i C, = 



F cU c m^ 

2pL^^b b 2c“ 

^ o o 



c d 

2^ a(u^t/c) 




(2.1S) 



in which is the reference velocity; U, the rate of deceleration of flow and z = f(^|^), 
i.e., the transformation given by Eq. (2.4). Equation (2.18) may also be deduced 
directly from the generalized Blasius equation. It is important to note that the force 
calculated from Eq. (2. IS) includes the effect of the rate of change of circulation 
between two successive time steps. Thus, it may be smaller or larger' (depending on the 
sign of F) than the force calculated through the integration of the instantaneous 
differential pressure Eq. (2.16). This is because of the fact that the instantaneous 
pressure depends only on the prevailing flow conditions and does not account for the 
rate of change of total circulation between successive time steps. In the calculations to 
follow and c are taken as unity for sake of simplicity. 

F. METHOD OF CALCULATION 

The methods used in the past in the determination of the vorticity flux from 
sharp-edged bodies may be roughly classified into two broad categories. The first of 
these involves the use of variable nascent vortex positions (see e.g., Sarpkaya 1968, 
1975) and the second, the use of fixed nascent vortex positions (see e.g., Clements 
1973-1975). 
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In the present analysis the method of variable nascent vortex positions is used. 
To e.xplain the computational details of the method let us consider a particular time t 
after the start of the motion and assume t to be sufiiciently large so that there are a 
number of vortices in the wake. Then the appearance and convection of the vortices 
proceed as follows: 

(1) Select a vortex position along the radial line defined by 0 = 117.72° (see 
Mostafa 1987). The very first location is taken r = 1.1 . 

(2) Calculate the strength of the nascent vortex which satisfies the Kutta 
condition. This is an exact solution and requires no iteration. 

(3) Place the nascent vortex at the corresponding points in the circle and physical 
planes and calculate the tip velocity. 

(4) Calculate a new nascent vortex strength from 0.5(V^-^ - V 2 ‘)At where is the 
tip velocity and V,. the average of three velocities along the radial line in the 
physical plane, i.e., at r = 1.05, 1.10 and 1.15. 

(5) Compare the newly calculated circulation with that obtained from the Kutta 
condition. If the difference between the two circulations is less than 0.001 
procede to the next step. If the said difference is larger than 0.001, carry out 
an iteration on the radial location of the nascent vortex as many times as 
necessary until the above condition is satisfied. If the circulation calculated 
from the Kutta condition is larger than that calculated from the tip velocities, 
the vortex must be moved towards the cylinder and vice versa. Also, each 
time the direction of the motion of the nascent vortex is changed (inward or 
outward), the marching distance is halved in order to accelerate the 
convergence of the two circulations. 

(6) Calculate the velocity induced at the center of all other vortices; 

(7) Convect the two nascent vortices with a velocity 0.5(Vj + V, ) for a time 
interval At (note that the vorticity is convected with the average velocity of 
the shear layer). 

(8) Convect all other vortices for the same time-interval At using a second order 
scheme given by 



z(t4-At) = z(t) -h 0.5 [3z(t) - z(t- At) ] At 



(2.19) 



in which z = u -f- iv. 

(9) Remove the vortices from the calculation whenever they come nearer than 
0.05 to the camber in the physical plane (e.xcept the first 20 vortices from the 
tip); 

(10) Coalesce the same sign vortices with a separation of less than 0.05 (in the 
physical plane, except the first 20 vortices); 
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(11) Calculate the tangential velocities and pressures on the inner and outer faces 
of the camber. Determine the drag and lift forces through the integration of 
pressure and through the use of the rate of change of impulse. Make plots of 
suitable variables (e.g., velocity distribution near the tip, variation of nascent 
vortex circulation with time, evolution of the wake, etc.); 

(12) Check the flow conditions to determine the state of the calculations: 

(a) IfVj - > 0.1 repeat the foregoing steps; 

(b) Stop the introduction of nascent vortices if 0 < Vj - V, < 0.1 and 
return to step No. 6; 

(c) If V 2 > Vj switch the angular positions of the nascent vortices to their 
image points. Calculate V., as the average of the three velocities, at the 
upstream side of the tip of the camber, at three radial locations (0.95, 
0.9, and 0.85) and repeat the foregoing steps; and 

(13) Make plots of the variations of various flow parameters (e.g., tip velocity, 
nascent vortex circulation, evolution of the wake, force coefficients, etc.) and 
terminate the run. 

The foregoing steps are quite general and can be used for any camber, provided 
that the optimum points of placement of the nascent vortices are determined through 
the use of a method developed by Mostafa (1987). 
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in. DISCUSSION OF RESULTS 



A. NUMERICAL AND PHYSICAL EXPERIMENTS 

The calculations were carried out for a time-dependent normalized velocity given 
by 



U t 

— =1 for T = -2— < 9.72 

*^0 



and 



— = 0.97T-0.05T^-3.70 
L' 



and 



(3.1) 



(3.2) 



A = 0.97-0.10T (3.3) 

in the interval 9.72 ^ T ^ 1 1.48 and 

U = -0.3423T + 0.0072T^ + 3.82 (3.4) 

and 

A = -0.3423 + 0.01445T (3.5) 

in the interval 1 1.48 ^ T. 17.95. For T larger than 18, the velocity and acceleration 
are zero. These velocities and accelerations correspond to that encountered in a series 
of experiments carried out in a vertical water tunnel. A detailed description of the 
equipment and procedures is given by Sarpkaya and Ihrig (1986) and will not be 
repeated here. Evidently, the calculations can be carried out for any specified variation 
of the velocity. For the case under consideration, the flow begins to decelerate at T* 



= 9.72 and the velocity of the ambient flow reduces to zero at about T* = IS, (see 
Fig. 3.1) 

The computer program provided, at times specified, the positions of all the 
vortices, the rate of shedding of vorticity from the tips of the camber, the velocity 
distribution on the upstream and downstream faces of the camber, the differential 
pressure distribution, and the force coefficients. 

Figures 3.2 through 3.5 show, at T* = U^t/c = 6.00, the evolution of the wake, 
the differential pressure distribution, and the velocities at the upstream and 
downstream faces of the camber. Figures 3.6 through 3.17 show similar plots of the 
wake, pressure, and velocity at larger times. These and other figures (not reproduced 
here for sake of brevity) show that the characteristics of the flow develop symmetrically 
prior to the onset of deceleration (T* < 9.72) and the differential pressure is positive 
everywhere (i.e., the pressure inside the camber is larger than that outside). 

Following the onset of deceleration (see e.g.. Figs. 3.8, 3.12, and 3.16), the 
differential pressure near the axis of the camber becomes increasingly negative. The 
significance of this result is that had the model been flexible (as in the case of a 
parachute) the central part of the camber would have collapsed as a result of the 
particular deceleration it is subjected to. 

For T* larger than about 13 (for the ambient flow under consideration), the 
velocities induced at the downstream edges of the camber by the large vortices moving 
sideways and towards the camber give rise to oppositely-signed vorticity. This, in turn, 
leads to the rapid growth of the secondary' vortices (see e.g.. Figs. 3.14 and 3.15 at T* 
= 21.275). Attention is drawn to the fact that the differential pressure shown in Fig. 
3.16 is negative over a large central portion of the camber. 

Figure 3. IS shows the variation of the circulation of the nascent vortex with T*. 
As expected, the vorticity flux is quite large at the start of the motion. It decreases 
gradually prior to the onset of deceleration and rapidly thereafter. For T* larger than 
about 12, i.e., when the secondary vortices come into existence, the circulation becomes 
negative. It reaches a minimum value when the ambient velocity and acceleration 
reduce to zero (at about T* = 18). The subsequent motion of the primary and 
secondary vortices increases the tip velocity and hence the strength of the nascent 
vortices. Nevertheless, the circulation remains negative, i.e., only the secondary' 
vortices continue to receive additional vorticity. 
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Figure 3.19 shows the velocities Vj and V, as a function of T*. The tip velocity 
Vj decreases from an initially large value of about 3 to a nearly constant value of 
about 1.4 just prior to the onset of deceleration. Subsequently, Vj decreases rapidly 
during the period of deceleration and prior to the inception of the secondary 
separation. Then increases to about 2 because of the backward motion of the large 
vortices near the tips of the camber. Finally, Vj decreases once again as the primaiy 
and secondary vortices move sideways and away from the tips of the camber due to 
their mutual induction (see Fig. 3.15). 

The variation of V, with T* is significant only during two, relatively short, time 
intervals; at the start of the motion and at the start of the deceleration. These are the 
periods during which the vorticity flux changes rapidly in order to maintain the Kutta 
condition. During the remainder of time V, is negligibly small, as expected on the 
basis of the pioneering experiments of Fage and Johansen (1928) with steady flow over 
various types of bluff bodies. 

Figure 3.20 shows the variation of the drag coefficient as a function of T*. It is 
calculated through the the integration of pressure about the camber. Figure 3.20 also 
shows that rises rapidly (due to the rapid accumulation of vorticity in the growing 
vortices) and begins to decrease as the vortices develop under the influence of a 
constant ambient velocity. Then the force decreases sharply at the onset of 
deceleration and goes through zero near the middle of the deceleration period (T* = 
14). The force acquires its largest negative value towards the end of the deceleration 
period. Subsequently, the force gradually decreases to zero. 

Figure 3.21 shows a comparison of the calculated and measured drag coefiicients. 
In general the agreement between the calculated and the measured drag coefficient is 
quite good. The differences are primarily due to the fact that the diffusion of vorticity 
has not been taken into account in the numerical analysis. It is possible to introduce a 
small artificial reduction in circulation in order to bring the calculated and measured 
values into closer agreement. This has been avoided in the present analysis in order to 
keep the discrete vortex analysis as pure and simple as possible. 

Finally, Fig. 3.22 shows a comparison of the normalized drag forces for various 
periods of the initial steady flow, prior to the onset of deceleration. It is clear that the 
drag forces beyond the point of deceleration are nearly identical. In other words, the 
force acting on the camber is not materially affected by the duration of the ambient 
steady fiow within the range of the parameters encountered in the present study. 
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Figure 3.1 Variations of the velocity and acceleration with T* 
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Figure 3.2 Position of vortices at T* = 6.00 
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.3 Instantaneous velocity of the vortices at T’’* = 6.00 
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Figure 3.4 Differential pressure distribution at T* = 6.00 
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Figure 3.5 Velocity distribution on both faces of the camber at T’* = 6.00 
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Figure 3.6 Position of vortices at T'*' = 12.275 
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Figure 3.8 DifTerential pressure distribution at T* = 12.275 
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Figure 3.9 Velocity distribution on both faces of the camber at T* = 12.275 
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Figure 3.10 Position of vortices at T'*' = 13.775 



35 



o 

cn 




LO 

O 



0 

~ CO 

1 

o 

cn 

I I ^ \ I I I \ I I I \ r I 

0‘9 O'S 0‘V 0‘C O'S 0‘I 0‘0 O'l- O'Z- 0‘£- O'V- 0'3- 0‘9- 

0/A 



Figure 3.11 Instantaneous velocity of the vortices at T* = 13.775 
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.12 DifTerential pressure distribution at T* = 13.775 
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Figure 3.13 Velocity distribution on both faces of the camber at T* = 13.775 
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Figure 3.14 
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Figure 3.15 Instantaneous velocity of the vortices at = 21.275 
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Figure 3.16 DiiTerential pressure distribution at T* = 21.275 
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.17 Velocity distribution on both faces of the camber at T* = 21.275 
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Figure 3.18 Variation of the circulation of the nascent vortex with T’” 
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Figure 3.19 Variation of the velocities and V, with T* 
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Figure 3.20 Calculated drag coefficient as a function of 
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Figure 3.21 Comparison of measured and calculated drag coefficients 
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22 Measured drag cceriicieats for various periods of initial steady flow 



47 



B. CONCLUDING REMARKS 

The results presented in this section have shown that the discrete vortex model 
can be used with confidence to predict the evolution of the wake about a cambered 
plate immersed in an arbitrary time-dependent flow. The drag coefficients resulting 
from the analysis and experiments agree reasonably well. This agreement can be 
improved with the introduction of a small circulation dissipation. 

The development of negative differential pressures near the central region of the 
camber is thought to be primarily responsible for the inception of the partial collapse 
of a parachute at high rates of deceleration. This phenomenon takes place even when 
the total drag force acting on the parachute is still positive. The sample analysis 
presented herein also shows that the negative differential pressure can cover a large 
region of the parachute and even result in negative drag. The basic idea emerging from 
the analysis reported herein is that the designs which incorporate into them the idea of 
delaying or preventing the return of the shed vortices to the canopy (e.g. .porosity 
management, change of deceleration history, parachute shape, dissipation and/or 
destruction of the organized w^ake) will be the ones which could avoid the collapse 
phenomenon. Extensive analysis and small scale experiments coupled with a few 
judiciously selected field tests may help to arrive at practically and phenomenologically 
sound parachute designs. 
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